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To complement hydrodynamic studies of the tidal disruption of the star by a 
massive black hole, we present the study of stellar luminosity and its variations, 
produced by the strong gravitational field of the black hole during a close en- 
q . counter. By simulating the relativistically moving star and its emitted light and 

£j , taking into account general relativistic effects on particle and light trajectories, 

c3 . our results show that the black hole's gravity alone induces apparent stellar lu- 

minosity variations on typical timescales of a few r 5 /c (=5 sec ) to a few 
><j : 100 i g /c (~ 10 min 1( ^6m ) ; where v g — Gm^/c 2 . We discern different cases with 

respect to the strength of tidal interaction and focus on two: a) a star encoun- 
tering a giant black hole traces space-time almost as a point particle, so that 
the apparent luminosity variations are dominated by clearly recognizable gen- 
eral relativistic effects and b) in a close encounter of a star with a black hole of 
similar size the stellar debris is spread about the black hole by processes where 
hydrodynamics plays an important role. We discuss limitations and results of 
our approach. 



1. Introduction 



Motivation for our work comes from the presence of massive black holes in galactic nuclei 
and from the possibility that such black holes accrete material from their surroundings. It 
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was estimated (Gurzadyan & Ozernoy 1981; Rees 1990; Magorrian & Tremaine 1999; Syer 
& Ulmer 1999) that central black holes may capture stars from inner galactic regions at the 
rate from 10~ 3 to 10 -7 per galaxy per year. Such events would be particularly interesting 
in Galactic centre, where the observed X-ray flare (Baganoff et al. 2001) and measured 
motion of stars, down to only 17 light hours from the centre (Schodel et al. 2002), provide 
strong evidence that the central concentration of about 3xl0 6 M is indeed a black hole. 
In recent years UV and X-ray flares have been observed in the nuclei of NGC 4552, NGC 
5905, RX J1242.6-1119, RX J1624.9+7554 and others, for which it was concluded that tidal 
disruption of a star by a massive black hole provides the best explanation (Renzini et al. 
(1995), Komossa & Bade (1999), Grupe et al. (1999), Gezari et al. (2003)). 

The interaction of a star with a black hole has been studied previously by other authors 
(Rees 1988; Carter & Luminet 1985; Luminet & Marck 1985) with a number of detailed 
hydrodynamics simulations (Laguna et al. 1993, Khokhlov et al. 1993a, 1993b, Kochanek 
1994, Fulbright et al. 1995, Marck et al. 1996, Diener et al. 1997, Loeb & Ulmer 1997, 
Ayal et al. 2000, Ivanov & Novikov 2001, Ivanov et al. 2003) with emphasis on stellar 
structure during the encounter with the black hole and longterm evolution of stellar debris. 
Nevertheless, none of these studied the luminosity variations occurring to the star in the 
vicinity of the black hole. In order to be complete, such study should include stellar hydro- 
dynamics in full general relativity, modeling of radiation processes in the disrupted star and 
relativistic effects on the emmitted light. Due to the complexity of the subject, we do not 
attempt to study all these effects in full here, but we wish to complement hydrodynamic 
studies by previously mentioned authors. Therefore we limit our attention in this paper to 
effects on star's luminosity induced solely by the gravity of the black hole, as we expect that 
relativistic effects alone might produce interesting luminosity phenomena. We simulate the 
disruption and the appearance of the star during close encounter as it would be seen by a 
distant observer and make a comparison of some results in our model with those obtained 
by hydrodynamic simulations. 

The model of the star used in our simulations depends on the expected strength of the 
tidal interactions between the star and black hole. Tidal disruption of the star with mass M* 
and radius R* occurs only if the star approaches the black hole to within its Roche radius: 




(1) 



which, expressed in units of the black hole's gravitational radius r g = Gm^jc 



• 2 , reads: 




(2) 
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where p Q and are the average densities of the Sun and the star. It is convenient to 
introduce the dimensionless Roche radius penetration factor f3 = r R /r p , where r p is the 
periastron distance of the star with respect to the black hole. The Roche penetration factor 
of a black hole grazing orbit is obviously: j3 gr = rn/(2r g + R*) = 1Zr/(1 + It is 

shown in the Appendix that this penetration factor crucially determines the strength of 
tidal interaction, i.e. the amount of work the tidal forces do on the star. We show (eq.A23) 
that tidal work can be approximated by: 

W ude ~ e\(3) = M*c 2 ^4 ^ ( 3 ) 

r p r p r p 

where e{(3) can be thought of as an effective eccentricity of the star at the periastron. If 
the Roche radius penetration factor is large, e may grow to values of order 1, bringing Wude 
to values comparable to a sizeable fraction of M*c 2 . Thus, the tidal interaction becomes 
overwhelmingly strong for large f3. Such an extreme scenario occurs for grazing interactions 
only if the size of the star is comparable to that of the black hole (see Appendix). We clasify 
grazing tidal interactions as follows: 

• m bh /M* <C 1: the Roche radius is smaller than the radius of the star, it follows that 
the Roche penetration factor is less than 1. As a consequence e 2 <C 1 and the star as 
a whole does not suffer large perturbations, even if the black hole pierces the star and 
accretes a small part of its mass along the way. 

• m bh /M* ~ 1: the Roche radius more or less equals the radius of the star (eq. 1) and, 
unless the star is very unusual i?* r g , so that the value of the Roche penetration 
factor Pg r ~ I. For such a f3 gr , e(/3 — 1) ~ 1 (c.f. Appendix A) and eq. 3 predicts 
that the tidal energy is of order lO^M^c 2 , which is a typical internal energy of a solar 
type star. Thus, the tidal energy is just about large enough to completely distort 
the star; the interaction may trigger violent hydrodynamic phenomena, possibly even 
a supernova. The most important phenomena governing the appearance of the star 
during such an encounter are hydrodynamic in nature, since the strong gravity region 
about the black hole has a much shorter range than is the size of the perturbed star. 
Hydrodynamics governs the appearance of the phenomenon and, therefore, such an 
event does not directly reflect general relativistic effects in strong gravity environment. 

• m bh /M* ~ (c/f e ) 2 , where v e = (2GM*/ R*) 2 is the escape velocity from the star: the 
black hole radius is comparable to the size of the star; if the star is not very unusual, its 
escape velocity is much less than c, so that according to eq. 1 the Roche radius is much 
greater than R* and consequently (5 gr ^> 1. In this case the tidal energy exceeds the 
internal energy by several orders of magnitude. A total and complete tidal disruption 
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takes place outside the black hole but in the region sufficiently close to the black hole 
for relativistic effects to play the mayor role in dynamics of the disruption. (Sect. 3). 

• (c/v e ) 2 < m bh /M* < (c/f e ) 3 : the black hole radius is larger that that of the star, 
but still smaller than Roche radius - (3 gr decreases with increasing mass of the black 
hole. The tidal energy before reaching the horizon is still comparable to the internal 
energy of the star. The release of tidal energy may well be sufficient to produce high 
energy shocks boosting stellar luminosity by many orders of magnitude. Yet, shocks 
moving with a few Mach are still much slower than the near speed of light the star is 
moving now. The star remains small with respect to the black hole along its way to 
the black hole. Such a stellar capture will thus very closely trace relativistic effects in 
the space-time, as it will be seen almost as a flashing up point particle on its way to 
doom. 

• mbh/M* > (c/f e ) 3 : the black hole is very much larger in size than the star (nibh > 10 8 M 
for a Solar type star), the Roche radius lies beyond the black hole's horizon, so it fol- 
lows (2) that the star is tidally disrupted only after crossing the horizon (t\r < 2r g ). 
Hence, the point particle approximation for the falling star is very good in the whole 
region outside the black hole. Since there is no agent to heat the star up, it is less 
likely for such an event to be noticed (Sect. 2). 

Here we discuss only the last three cases, since we find them interesting as a tool to 
study the strong gravity regions in the universe, as well as in view of supermassive black 
holes in galactic nuclei. 

2. Star encountering a m bh > 10 s M Q black hole 

We expect that a capture of a star by a giant black hole would most likely occur when 
a star in the cluster surrounding the giant black hole would be perturbed to a low angular 
velocity orbit with respect to the black hole. Therefore, such encounters will, most likely 
occur with the velocity characteristic for parabolic infall. During such an infall the star can 
not be significantly disrupted while outside the horizon, so with respect to a much larger 
black hole it can be treated as a point source of light whose appearance with respect to the 
far observer will be modulated by the Doppler shift, aberration bending and gravitational 
redshift. Two numerical codes were developed to calculate the apparent luminosity changes 
of the source falling in both, Schwarzschild type and Kerr type black holes. During the 
encounter of the star with a giant black hole the star is simulated as a point source emitting 
monochromatic light of frequency z/ and intensity L , both constant in the frame comoving 
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with the source. As the source is moving along a parabolic orbit with a given orbital angular 
momentum, we trace light rays from subsequent points of the source's trajectory (separated 
by At=l r g /c in coordinate time) to the distant observers and calculate the apparent lumi- 
nosity with respect to them as a function of time. We would like to note, that these results 
are directly applicable also to luminosity and spectrum changes produced by orbiting blobs 
of material in accretion discs around black holes. 

Results for both types of black holes show (Gomboc et al. 1999) two characteristic 
timescales of luminosity changes, both determined by the gravity of the black hole. The 
first one displays the basic quasiperiod in luminosity and redshift changes as the star spi- 
rals toward the black hole. The quasiperiod very closely matches the orbital period of the 
source at the innermost stable orbit (50 r g /c for Schwarzschild black hole). The number of 
quasiperiods observed depends on the fine tuning of the angular momentum to the critical 
value. In Schwarzschild case the critical angular momentum is l/M*r g c = I =4 and the 
number of quasiperiods can be approximated as N p = 0.5 — 0.5Log(4 — /) for 3.9< I <4. 
The quasiperiods in Kerr case differ for prograde and retrograde orbits: for a maximal Kerr 
black hole (with rotation parameter a=0.998 r g ) and a star on a prograde orbit with angular 
momentum close to critical l + = 2(1 + — f-), the quasiperiod is ~ 13 r g /c, while for a 

star on a retrograde orbit with / close to critical l~— —2(1 + ^Jl + f-), the quasiperiod is m 
80 r g /c, both consistent with orbital periods at critical radii. 

The second time scale is considerably faster (of order lr g /c) and belongs to the rate 
of change of relativistic beaming direction with respect to the observer. For the black 
hole with mass iribh = 10 8 M Q , the corresponding timescales are ~ 10 hours and ~ 10 min, 
and for extreme nibh = 10 10 M Q , this time intervals are ~ months and ~ 10 hours. Since 
the luminosity and spectrum changes are caused by relativistic beaming and gravitational 
lensing, they are most evident to observers in the orbital plane of the star. The observers 
perpendicular to this plane see the source as slowly fading and then, as the source approaches 
the horizon, suddenly disappearing on a timescale of the order of ~ 10 r g /c. Comparing 
results for Schwarzschild and Kerr black holes, we find that luminosity curves (Fig. 1) are 
qualitatively similar, but timescales generally shorten for Kerr prograde orbits and become 
longer for Kerr retrograde orbits. Results show that within 5° of the orbital plane one may 
expect luminosity rise by a factor of a few 10, while the maximum Doppler plus redshift 
factor (i/obs/vo) is 1.8 for the Schwarzschild case and 2.2 for the maximal Kerr black hole 
case. 
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3. Star encountering a m b h ~ 10 5 - 10 6 M black hole 

3.1. Approximations, model and comparison with hydrodynamic results 

The capture of a star by a black hole of comparable size is a phenomenon where black 
hole's gravity plays the dominating role both on propagation of light as well as on propagation 
of matter belonging to the star. This property of the phenomenon is forcefully stressed by 
the fact that the tidal energy is many orders of magnitude larger than its gravitational 
binding energy and becomes a sizeable fraction of M*c 2 (eq. 3). Therefore, we build our 
approach on the work of Luminet & Marck (1985), who showed that in the vicinity of the 
black hole "particles of the star undergo a phase of approximate free fall in the external 
gravitational field, since the tidal contribution grows much larger than pressure and self 
gravitating terms". Therefore, we use a simple model, whereby the star is considered as 
undisturbed by the black hole (i.e. spherically symmetric), until it reaches the Roche radius. 
After crossing it, the black hole's gravity takes over and the self gravity and internal pressure 
are completely switched off. 

Further, we neglect hydrodynamic effects. This approximation is justified if the proper 
time elapsed between the Roche radius crossing and total disruption is short compared to 
the dynamic time scale of the star. For an estimate of the two timescales we take 

tr ~ ^r g -nl = (6nGg^, (5) 

where tr is estimated as the proper time elapsed during a radial parabolic infall from the 
Roche radius to the horizon 1 . Specifically, for a solar type star we obtain tr ~ 13 min, which 
is an order of magnitude less than the dynamic time scale r d m 3 hours. The ratio of the two 
times indicates that the amount of energy exchanged may not be quite negligible, but is small 
enough that it may be neglected in the first approximation. Further justification for such an 
approximation comes from results of hydrodynamic evolution calculated by Kochanek (1994) 
and Laguna et al. (1993). Laguna et al noted that "the qualitative features of the debris - 
including its crescent-like shape - can be reproduced by neglecting hydrodynamic interactions 
and self-gravity of the star", since the formation of the crescent is due to "geodesic motion 
of the fluid elements of the star in a Schwarzschild space-time which includes relativistic- 
induced precession of the orbit about the black hole" . This confirms previously mentioned 



1 Of course, tr is denned only for tr > 2r g , when tidal disruption takes place outside the horizon of the 
black hole. For nonzero angular momentum orbits tr is slightly, but not crucially longer. 
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findings by Luminet & Marck (1985), that black hole's gravity dominates in close encounters. 
Therefore, we argue that by neglecting hydrodynamic effects, we obtain in close encounters 
approximately the correct shape of stellar debris. 

Hence, our numerical model starts with a spherically symmetric star of radius i?* and 
mass M* consisting of N equally massive constituents (m* = M*/iV, N ~ 10 6 ) distributed 
randomly but in such a way that in the average their density distribution follows that of a 
star, which is approximated by the polytrope model with n=1.5. All constituents of the star 
start with the (same) velocity, corresponding to parabolic velocity of the stellar center of 
mass, which is placed at a distance TZr from the black hole. Subsequently the positions of 
free falling stellar constituents are calculated at later discrete times (ti) according to general 
relativistic equations of motion. 

To test the errors induced by these approximations, we tudy encounters of a M Q , R Q 
star with a 10 6 M Q black hole and compare our results on central density in the star (average 
density inside 0.01 i?*) with those obtained by Laguna et al (1993), Fulbright et al. (1995), 
Khokhlov et al (1993b), Ivanov & Novikov (2001). Fig. 2 shows the central density as a 
function of time with respect to the time of periastron passage as obtained by our model and 
by hydrodynamic simulations. The qualitative agreement between these results justifies the 
neglect of internal pressure in calculating the dynamics of disruption. The major difference 
seems to be in the precise timing of tidal compression: in our model the strongest compression 
occurs very close to periastron, in agreement with the results of Luminet & Marck (1985), 
while in most hydrodynamics simulations the central density peaks approximately 15-20 r g /c 
after the periastron passage. 

Our results on the shape of stellar debris during the close encounter also agree with 
results of Laguna et al. (1993), although at later times our crescent becomes considerably 
longer. 

3.2. On the luminosity of the star during the tidal disruption 

We consider the tidal disruption to be the phenomenon, where the work done on the star 
by tidal forces is comparable or greater than its initial internal energy. The tidal disruption 
is thus a violent nonstationary process that takes place in the vicinity of the black hole 
on a time scale that is considerably shorter than the stellar dynamic time scale (measured 
in proper time of the falling star). As the star is deformed into a long thread, the giant 
tidal wave deposits great amounts of energy which soon pushes gases into an out moving 
shock wave more or less perpendicular to the threadlike axis of the star. Thus, during the 
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disruption process several mechanisms play an important role: shocks, adiabatic expansion 
and cooling of disrupted material, possible explosions due to tidal squeezings as predicted by 
Carter & Luminet (1982, 1985), radiation driven expansion etc. These effects have no doubt 
important influence on the cooling and luminosity of the disrupted star, but we wish to 
stress, that as shown by Luminet & Marck (1985) gravity in general overwhelms other forces 
during the close encounter. So, since gravity of the black hole swings the star around on a 
time scale that is much shorter than any other time scale that may play a role, we believe, 
that as a first step to estimate the luminosity variations of the tidally disrupted star, we may 
use a simple model, which must in the first place correctly take into account the effects of 
dominating strong gravitational field of the black hole. As the disruption progresses and the 
hot stellar inner layers become exposed both by gravity and by shock waves, the luminosity is 
bound to rise due to higher effective temperature and due to higher effective area seen. The 
overall rise in luminosity depends on other partially competing mechanisms involved: while 
the expansion and cooling would tend to reduce it, it must nevertheless rise dramatically due 
to enormous work being done by tidal forces which drive shock heating and supernova-like 
explosions. The precise role of these mechanisms and their influence on stellar structure and 
evolution need detailed analysis, but is beyond the scope of this paper. 

Here we wish to make a step towards the complete solution by including in full only 
the most important ingredient defining the shortest time scales: the effects of black hole's 
gravity on the apparent variability of stellar luminosity. The standard stellar atmosphere 
model (Bowers and Deeming 1984, Carroll and Ostlie 1996, Swihart 1971) is not applicable 
in calculating the effective temperature of any surface element since, because of the highly 
dynamic structure, the fine details of atmospheric density, temperature and pressure profiles 
are not available, even more we can not predict in advance which part of the star is going 
at some future time belong to the atmosphere. So we are forced to apply a Monte Carlo 
model throughout the star by which the unperturbed star is modeled as a spherical cloud 
consiting of a large number (N) of identical constituents distributed randomly, but in such 
a way that their average density follows that of an n = 1.5 polytropic model (cf. section 
3.1). The constituents are optically thick and have an assigned temperature according to 
their position in the cloud, which again follows the temperature profile of the n — 1.5 
polytropic stellar model. The model "photospheric" temperature and model "luminosity" 
are calculated as the sum of spectral contributions from those cloud constituents that are 
seen by the observer, i.e. by those that are not obscured by constituents in overlaying layers. 
For the purpose of obscuration all the constituents are considered to have the same cross 
section a, so that a = 4ttRI/N', where the parameter N' is the number of constituents 
belonging to "the atmosphere" of the star. It is clear that, since for statistical reasons N' 
must be at least a few ten, and iV is limited by the computer power to a few million, the 
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ratio N' /N is much greater than the ratio Matm/M* in a real star. One could argue that 
the atmosphere could be made less massive by representing it with a larger number of less 
massive constituents. However, in the case of total tidal disruption the interior is mixed into 
the atmosphere during the late stages of disruption and the so introduced uneven opacity of 
stellar constituents would further complicate the interpretation of results. Thus we can not 
afford to make models with sufficiently opaque atmospheres and, as a result, our initial model 
" photospheric" temperatures are too high. We note, however, that the model photospheric 
depth is a function of N'/N, so by changing N, we probe the stellar atmosphere to different 
depths. In such a way an extrapolation to realistic opacities is possible. The consistency of 
such an extrapolation is checked on the initial spherically symmetric stellar model, where 
the Monte Carlo results can be directly compared with the theoretical atmospheric model. 
An example of such a comparison is shown in Fig. 3. It is clear that the depth of our 
model "photospheres" is some orders of magnitude too high, yet it is possible to extrapolate 
model "photospheres" to depths of realistic stellar atmospheres, since the temperature is 
a monotonic smooth function of depth. For evolved stages of tidal disruption there is no 
underlying theoretical model, so that we rely on extrapolated results of the Monte Carlo 
model. 

As the star moves along the orbit, images of the star with respect to the far observer are 
formed as follows: Photon trajectories and the time of flight between each stellar constituent 
and the observer are calculated with technique described in Cadez et al. (2003), Gomboc 
(2001), Brajnik (1999) and Cadez & Gomboc (1996). Only two trajectories connecting two 
space points are considered - the shortest one and the one passing the black hole on the other 
side, while those winding around the black hole by more than 2ir are neglected. (It has been 
shown before (Cadez & Gomboc 1996), that light following trajectories with higher winding 
numbers contributes less and less to the apparent luminosity.) The beam contributions are 
sorted into pixels with an area corresponding to the size of a, and tagged according to the 
arrival time. The intensity corresponding to a given pixel is then defined as the intensity 
corresponding to the ray with shortest travel time. Since light from deeper layers takes longer 
to reach the observer, this takes care of the obscuration of deep layers. The intensity of a 
contributing beam is calculated assuming that the corresponding stellar constituent emits 
in its own rest frame as a black body at its temperature. The apparent luminosity and 
effective temperature of the star as a function of time (with respect to the chosen observer) 
are calculated and successive stellar images, formed in this way, are pasted in a movie 2 . 

We divide our model in three parts: first we estimate the relativistic effects alone by 



2 Movies can be obtained at www.fmf.uni-lj.si/"' gomboc 
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simulating the luminosity variations of an iso-thermal star (i.e. star with T(r,t)=const.). In 
the next step, we consider the star with polytrope n=1.5 temperature profile and estimate 
the luminosity variations due to exposure of inner hot regions of the star. We first consider 
a simple case, in which the temperature of all stellar constituents is constant in time (no 
cooling or heating), and afterwards add a rough estimation of the effect of cooling of exposed 
stellar parts on the stellar luminosity. 

3.2.1. Effects of black hole's gravity 

To isolate the effect of gravity, we first compute luminosity variations of an iso-thermal 
star. The ensuing luminosity variations can be ascribed to: Doppler boosting and aberation 
of light, gravitational lensing and redshift (similar as for a point-like source in Sect. 2) and (in 
addition) the elongation of the star due to relativistic precession and due to tidal squeezing. 
Fig. 4 shows the obtained luminosity variations as a function of time for encounters with 
1=0 (radial infall), 1=4 (critical), 1=5 (r p =10 r g ) and 1=7 (r p =22.3 r g ) as seen perpendicular 
to and in the orbital plane. 

Results show that the maximal rise in luminosity occurs in the case of the critical 
encounter (1—4), where the overall luminosity rise due to elongation of the star is of about 
a factor of 20 (as seen by the observer perpendicular to the orbital plane, Fig. 4 above), 
while gravitational lensing and Doppler boosting enhance it up to about 40 times the initial 
luminosity (Fig. 4 below). Observers close to the orbital plane see most extreme variations: 
dimming of the receding star, its rebrightening as it emerges from behind the black hole and 
variations on short timescales of about 10 r g /c, which are due to lensing effects. Since the 
star and the black hole are comparable in size, the probability that they are aligned with 
respect to the observer, is high. When lensing takes place the relevant part of the star is 
imaged into an Einstein disk and the apparent luminosity increases manifold (Fig. 4c). 

3.2.2. Constant temperature debris 

Next, we consider the star with n=1.5 polytrope temperature profile and we assume 
that the temperature of stellar debris does not change with time. The model is obviously 
much too crude to rely upon its results regarding the spectral characteristics or even the 
absolute value of the emitted luminosity. The crude argument why this model may bear 
some resemblance to the true light curve is that shock wave released by the unbalance of 
gravity carries internal energy to the surface in such a way that the energy influx from the 
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interior temporarily compensates the radiation loss. 

The simulation shows that, as the inner hot layers of the star are exposed during the 
disruption, they contribute to the substantial rise in stellar luminosity, depending on the 
orbital angular momentum of the star (Fig. 5). The star on a low angular momentum orbit 
is completely captured by the black hole and produces only a short (~ 1 - 10 r g /c) flare before 
disappearing behind the horizon. On the other hand, the star with high angular momentum 
experiences only a slight distortion during the distant flyby with a resulting temporary (~ 
10 - 100 r g /c) slight increase in luminosity. 

The most dramatic is the encounter of the star on the critical angular momentum orbit 
(7=4), during which half of stellar constituents are swallowed by the black hole and the other 
half escapes. During this process, the star is totally tidally disrupted in such a way that the 
higher angular momentum material rapidly lags behind the stellar debris with lower angular 
momentum, which produces a long thin spiral (Fig. 5.). Outer layers of the star are stripped 
off in a time of the order of 100 r g /c, the depth to the hot inner core decreasing together 
with self gravity. In our crude model this is seen as decreasing optical thickness and the 
exposure of the hot inner core; the luminosity rises steeply. The spectrum of the debris is 
dominated by the emission of the innermost exposed layers and as long as shock waves are 
building up, i.e. until cooling sets in, these lead to X-rays. 

Some luminosity peaks arise from the effect of tidal compression in the direction per- 
pendicular to the orbital plane of the star, which in our model for a short time exposes the 
interior of the star. Such peaks are evident in Fig. 5 above c and c', and these two com- 
pressions are in agreement with multiple tidal squeezings predicted by Luminet & Marck 
(1985) and confirmed by Laguna et al. (1993). In our model they produce luminosity peaks 
lasting about ~ 5 r g /c. As mentioned earlier, Carter & Luminet (1982, 1985) predict, that 
thermonuclear explosion may occur at this moment. 

The scale of the luminosity rise in Fig. 5 is rather uncertain due to neglect of hydrody- 
namic effects 3 and also due to our poor atmosperic model (section 3.2.). For the critical tidal 
disruption of the Sun the extrapolation of our model would suggest the total luminosity to 
rise to about 1O 13 L (mostly in X-rays), which accentuates the extent of tidal disruption, 
but also sends a warning that by that time our constant internal energy model assumption 
ceases to be valid. As suggested in section 3.2., we calculated a range of models with N 
between 10 3 and 10 6 and extrapolated the results to realistic atmosperic depths. These nu- 



3 For simplicity, we assume that all the tidal energy is transformed into the kinetic energy of the tidal 
wave; the portion of kinetic energy that may go into heat is neglected, therefore, we expect that the actual 
available luminous energy during such a tidal disruption may be higher than the one given by our models. 
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merical results suggest that, at least for the critical disruption, the average temperature and 
size of the final crescent to which the star is deformed is roughly independent of N. Thus 
we tested the idea that tidal disruption exposes or mixes up by shearing the envelope of the 
star to a certain depth R c , which we define as the depth in the undisturbed star, down to 
which the average T 4 is equal to the average T 4 of the final crescent. In this way we estimate 
(independent of N) that for critical 1—4 and n =1.5, R c is about 0.25 i?*, while for n = 5 
we get R c =0.1 i?*. For a close ffyby with 1=5, R c is about 0.7 i?* and 0.5 /?* for n =1.5 
and n=5 respectively. We may also, as an example, estimate the luminosity of a Solar type 
star during the bright critical stage of total disruption on a 1O 6 M black hole as follows: the 
steep luminosty rise (c.f. Fig. 5) has a time scale between 30r g /c to 100r 3 /c, which is about 
2.5 to 8 minutes. Assuming that the initial thermal energy contained in the exposed layers 
(~10 48 erg) is radiated away on this time scale, the critical luminosity would be of the order 
5 to 15xlO n L . 

After the debris is spread and starts moving away from the black hole, the physics of 
tidal disruption is no longer dominated by black hole's gravity. The physical conditions 
in stellar debris, the physics of radiation processes, hydrodynamics etc. take over and the 
ensuing processes go beyond the simulation presented here. 

3.2.3. Cooling of stellar debris 

In general, the temperature inside the star may change due to various mechanisms, 
already mentioned. To get an idea of how they might affect the light curve, we again model 
the cooling in two very approximative ways: 

(a) exponential cooling of exposed stellar layers with different characteristic times: r=l 
r g /c and 10 r g /c, 

(b) cooling of exposed stellar layers as due to their own black body radiation in the 4ir 
solid angle. 

Results presented in Fig. 6 show, that if the cooling were very efficient - with timescales 
of 1 r g /c, the luminosity rise would be quite short and modest. 

4. Conclusion 

Stellar encounter with a massive black hole can be a very energetic event, with energy 
released and luminosity variations depending primarily on the relative size of the star com- 
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pared to the black hole. We note that the tidal interaction energy may rise to as high as 10% 
of the total mass-energy of the captured star, which is available when the star is comparable 
in size to the size of the black hole. This size ratio is also critical as to the nature of the 
disruption. 

In this work we focused on gravitational phenomena and showed that: 

• a) A critical capture of a "pointlike star" is characterized by a series of quasiperiodic 
apparent luminosity peaks with the quasiperiod 50 r g /c for a Schwarzschild black hole 
and 13 and 80 r g /c for an extreme Kerr co- and counter-rotating case respectively (Fig. 
1). This translates into 6.9 hours x 1( ^ , 1.8 hours x 1( ^ and 11.1 hours x , 
respectively. If a "pointlike star" would be a planet falling to a 3.6 x 10 6 M Q black 
hole in Galactic centre respective quasiperiods would be 15 minutes, 3.9 minutes and 
24 minutes. 

• b) The sharpness, the amplitude of quasiperiodic peaks and the amplitude of the 
Doppler factor is more pronounced for observers in the orbital plane as compared to 
those perpendicular to this plane. The highest value for the Doppler factor is 1.8 for 
the Schwarzschild and 2.2 for the extreme Kerr black hole. 

• c) The number of quasiperiodic peaks (N p ) depends on the closeness of the orbital 
angular momentum (/) to the critical value 1=4 and can be approximated as N p = 
0.5-0.5Lo#(4-/~). 

• d) An extended star may be approximated as a collection of point particles when 
heading toward the complete tidal disruption. The shape and the density of the debris 
calculated in this approximation compares well with more sophisticated hydrodynamic 
calculations (cf. Section 3). 

• e) Model light curves for critical tidal disruption of a star of the same size as that of the 
black hole (Fig. 4, 5, 6) calculated for different heuristic models show similar temporal 
characteristics which display very rapid (on time scale of order 10 r g /c) luminosity 
variations by a few or even many orders of magnitude, while the quasiperiodicity is 
no longer pronounced in such a process. Light curves describing a critical capture are 
very rough and can not be momentarily calibrated in flux. They are presented as they 
produce the extremely short time scale phenomena characteristic of the strength of 
black hole's gravitational field, which will persist in the future more elaborate models 
of tidal disruption. 
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A. The virial theorem and tidal energy 

In order to estimate the amount of heat and kinetic energy deposited to the star by the 
tidal wave, it is useful to follow the steps of the derivation of the virial theorem. Consider the 
some 10 60 nuclei and electrons making up the star as representative point particles making 
up the ideal gas of the star. Each of the particles with mass m ; (i = 1 . . . 10 60 ) moves 
according to Newton's law (We will folow the more transparent classical derivation, which 
is sufficient for order of magnitude arguments.): 

3& 



The black hole has been placed at the origin from where the position vectors are reconed. 
F£ models the force taking place during particle collisions. It obeys (the strong version of) 
the third Newton's law, and since in the ideal gas approximation collisional forces act only 
at a "point", the energy connected with the potential of these forces can be neglected. The 
second term on the right describes the gravitational interaction among the constituents of 
the star and the last term represents the gravitational force of the black hole. It is convenient 

to define the center of mass position vector R = rriif^j /M*, so that r*j = R + r* and 

^mjf* = 0. Summing equations Al over all i, one obtains the center of mass equation of 
motion in the form: 

M J = -G^^R-5Gm bh ^^R + 2Gm bh ^ + 0(l/R 5 ), (A2) 

where Q is the quadrupole moment tensor of the mass distribution with respect to the center 
of mass of the star defined in the usual way as: 

i 

Terms of 0(1/ R 5 ) and higher will henceforth be neglected. If the star is deformed in a 
prolate ellipsoid with the long axis in the direction n, Q can be written in the form 

Q = 3qhh-ql , (A4) 

with q beeing positive and proportional to the eccentricity of the ellipsoid. Here fjfj stands 
for the diadic product of the respective vectors and I is the identity matrix. 

The angular momentum of the star (/), which is a conserved quantity, can be split into 

the orbital (l Q = MR x R) and spin part (l s = X)z m ^i x ^ ne time derivative of the 
orbital part follows from eq.A2 and when A4 applies, it can be written as: 

t = 6G^0-(Rxn)(R.n) (A5) 
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The sum of scalar products of equations Al by r*j gives the energy conservation law. We 

1 u2 

split the kinetic energy of the star into the center of mass part ^M^R and the internal kinetic 
energy part Wi nt = ^ i m i^i ■ Using eq. A2 and neglecting the collisional interaction 
energy, we obtain the conserved energy E in the following form: 

ei 1 d n m bhM* m bh R.Q.R 

E = -M*R - G — G — 5 + W int + W G , (A6) 

where Wq is the self gravitational energy of the star (Wq — —\ J2i G ^*™!, ) 

J i I j i I 

Finally, we obtain the equivalent of the virial theorem by we multiplying eqs.Al by r* 
and summing over all %. The result can be rearanged into the transparent form : 

W mt + \W G = - G 1 ^^ + \j , (A7) 

where J = Yli m i r 'i ' ■ For a star in hidrostatic equilibrium, the right-hand side vanishes and 
the total energy of the star W tot = Wi nt + Wq = — Wi nt . If the star is not in hidrostatic 
equilibrium, the right hand side of eq.A7 can be considered as the energy imbalance - if 
it is more than Wi nt , it is sufficient to completely disrupt the star on a time scale r<j. An 
exact evaluation of this energy imbalance is beyond reach in this simple analysis, however, 
a simplified model offers some clues. 

Consider an idealized case of an "incompressible star" flying about a massive black 
hole. From the point of view of the star, gravity is exerting a tidal force squeezing it in 
the plane defined by the temporary radius vector and the orbital angular momentum and 
elongating it perpendicular to this plane. The tidal force acts to accelerate the surface of 
the star with respect to the center of mass, but it must also act against rising pressure 
and internal gravity. Thus, roughly speaking, the tidal force does work in pumping kinetic 
energy into the tidal wave, but also in loading the gravitational potential energy which acts 
as the spring energy driving oscillation modes of the star. Consider small tidal distortions. 
In this case quadrupole deformations are dominant, so that the deformation field (U) of the 
incompressible star can be described as a linear combination of 5 degenerate quadrupole 
modes: 

5 

U = ^a fe II fc . (A8) 

k=l 

Here are modal base vector fields that can be expressed as gradients of quadratic poly- 
nomials in coordinates x', y', z', obtained by multiplying spherical functions Y 2m (9',(f)') by 



4 Note that W int comprises both the kinetic energy of thermal motion and the kinetic energy of bulk 
motion in the tidal wave. 
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r' 2 and identifying x' = r' sin 6' cos 0' etc, and a,k(t) are modal amplitudes. In the coordinate 
system where the z' axis is normal to the orbital plane and x' points from the periastron to 
the black hole, only three amplitudes are excited and the corresponding modal base fields 
are: 

n 5 = -y — (-x,-y,2z) n 1 = -y— (x,-y,0) II 2 = — ( y, x, 0) (A9) 

These deformations lead to the following quadrupole moments: 

/ ai 2 + a 2 2 — 2 ai — a^, 2 , — 2 a 2 05, \ 
Q = — -2 v / 3a 2 a5,ai 2 + a 2 2 + 2v^a 1 a 5 -a5 2 ,0 (A10) 
V 0,0,-2( ai 2 + a 2 2 -a 5 2 ) / 

As long as tidal modes can be considered roughly independent, their dynamics can be derived 
from the Lagrange function L = T — U with the kinetic energy (T): 

5 5 p q 5 

T = / mdiUkUidV = ^M*RlJ2 d * (All) 

fc=i i=i i=i 

and the potential energy ([/ - the deviation od self gravity from the equilibrium value in 
undeformed state): 

U = —M^Rl^W , (A12) 
i=i 

where u) q is the resonant frequency of quadrupole modes. For a star consisting of a self 
gravitating incompressible fluid we obtain 

64 

u? q = -GMJRI (A13) 
Generalized forces exciting these modes are (Goldstein 1981): 

F fc = JpU k G^.(I-3^).PdV' , (A14) 

Let us calculate these forces in the specific case when one can assume that R(t) represents 
a parabolic orbit. We express the compontents of R as 

R(t) = i?(t){cos^(t),sin^(t),0} , (A15) 

where 

R(t) = r p /sin 2 -^(t) (A16) 
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and ip{t) is the true anomaly obeying the Kepler equation: 



^T* = -cot^(l + icot 2 ^) . (A17) 

With this, and using A9, the integrals in A14 can be evaluated to obtain the nonvanishing 
generalized forces: 




m bh I 3 „ 2 • 6 1 



COS 2ip 

" 3G l^V^ MXsin °2^l Sin2 ^ 1 (A18) 
Finally we write down the Euler-Lagrange equations of motion (^J^ — = Qk) for modal 

/ 2r 3 

amplitudes. After introducing the characteristic time tf — y G J/ bh and the dimensionless 
time r = t/tf, they can be cast into the dimensionless form: 

^ + (u q t f ) 2 ai = h(r) , (A19) 
where the dimensionless forces fiij) are functions of ipij) only: 

' 71 sin 6 -tj} sin 2ip I (A20) 





4 V 5 2 



Thus, the only trace of parameters of the tidaly interacting system is left in the factor u q tf, 
which is 27r times the ratio of the characteristic fly-by time about the black hole and the 
period of quadrupole modes. It is useful to note that, using eq.l and A13, this product can 
be written as 

u q t f = 8^2/5(r p /r R )l = 8^2/5(1/(3)1, (A21) 

i.e., it is inversely proportional to the power of Roche radius penetration depth. In the case 
of a distant fly-by uo q tj 3> 1, so it follows from eq.(A19) that a, L = fi(oj q tf)~ 2 oc y^, which 
is the familiar result often used with Earth tides. Note, however, that for deep penetrations 
of the Roche radius u q tf < 1, and thus the (dimensionless) generalized forces fi(t) become 
large at frequencies that are resonant with u q . 

We calculate the total work done by tidal forces on the system of normal modes during 
the whole fly-by process by noting that it can formally be expressed as the change of the 
Hamiltonian H(t) = T + U during the process (neglecting damping of normal modes). 
Initially the quadrupole system starts in the undisrupted state with H(t — > — oo) = 0, and 
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ends in a state of excited quadrupole modes with Wud = H(t — > oo) 5 (i.e. for t >> tf): 

q 5 5 »oo 

W m e = — M^Vlim^ + ^a, 2 ) = V / F^M* (A22) 
Solving equations A19 with the retarded Green's function, this can be written in the form: 

W ude = iG^^^U^tf)? (A23) 
p i=i 

where 

1 r°° 

m) = 7^ J J {T)einrdT (A24) 

We note that Wude can be written in the form Gm bh q/rp, where q = M^R^e 2 and 
according to (A23) 

e 2 = \ £|/;K*/)| 2 (A25) 
i=i 

can be thought of as an effective eccentricity of the star at the periastron. Fig. 7 shows that 
e can reach values of order 1 if a fly-by is comparable to the dynamic time-scale of the star. 
Note however that for deep Roche radius penetrations our first order perturbation model no 
longer applies; closer analysis shows that the model is aplicable for u q tf > 1 i.e. for (3 < 3 
(eq. A21) 6 . 

Now we are in the position to estimate the high value of the right hand side of eq.A7 for 
this simple parabolic infall of an incompressible star. The left hand side starts at zero, when 
the star is still far from the black hole. As time goes on, the internal kinetic and potential 
energy change as the energy of tidal modes, so that the left hand side is greatest when all 
the tidal energy is in the kinetic energy of the wave. Thus, the maximum value, which is 
also the maximum value of the right hand side equals W^de- 

Even if the above analysis is valid, strictly speaking, for an incompressible star and 
in the approximation of independent (small amplitude) tidal modes, it does suggest the 
qualitative conclusion that the tidal interaction depends crucially on the ratio period of the 
fundamental mode versus typical fly-by time (ou q tf) and does become resonant if the fly-by 



5 Assuming the tidal kick did not break up the star by imparting higher than escape velocity to surface 
layers. 

6 We note that for 1 < (3 < 3 the tidal energy is proportional to f3 2 since e 2 oc 1/(3. This is in agreement 
with result of Lacy et al 1982 and Carter & Luminet 1983. 
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time is less than the period of the fundamental mode. The energy deposited into the star by 
the tidal interaction can be of the order G mbhM 3 * R * = M*c 2 r g Rl/rl, which may surpass the 
absolute value of the internal gravitational energy of the star by many orders of magnitude 
if r p , and r g happen to be of the same order. 

We thank the anonymous referee for constructive criticism which helped us improve the 
text. We acknowledge the financial support of the Slovenian Ministry of Science, Educa- 
tion and Sport. AG also acknowledges the receipt of the Marie Curie Fellowship from the 
European Commission. 
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Fig. 1. — Luminosity and frequency shift during the infall of a solar type star into a giant 
black hole nibh > 10 8 M Q . a) infall with orbital angular momentum I close to critical 1=4 
into the Schwarzschild black hole, as observed perpendicular to the orbital plane; b) same 
event observed in the orbital plane; c) infall of the star on prograde orbit with / close to 
critical l + into the Kerr black hole, as observed in the orbital plane; d) infall of the star 
on retrograde orbit with I close to critical l~ into the Kerr black hole, as observed in 
the orbital plane. The color code in frequency diagram corresponds to spectrum intensity 
(in units of the initial intensity of the primary image). (Full resolution images available at 
www.fmf.uni-lj.si/~ gomboc) 
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lime[r e /c] time [r B /c] [imc Ir./c 

Fig. 2. — The central density in the star as a function of time during close encounter, for 
polytrope n=1.5: (3=5 (1=5), (3=10 ([=4.08) and n=2, r] = 0.1 ( rj = (^) 1/2 (^) 3/2 )- Solid 
curves are from our simulations, dashed are from Laguna et al. (1993), dotted from Fulbright 
et al (1995), dot-dashed from Ivanov & Novikov (2001) and short-dashed from Khokhlov et 
al. (1993). Time is measured from the periastron passage. (Full resolution images available 
at www.fmf.uni-lj.si/~ gomboc) 
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Fig. 3. — Left: The model atmosphere depth in the initial (spherically symmetric) star: 
polytrope temperature profile (line) and " photospheric" temperatures and corresponding 
depths for different N in our model (symbols). Right: "Photospheric" temperature as a 
function of time during total tidal disruption given by our model for different N. 
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Fig. 4. — Isothermal star with R Jf =2r g during the encounter with critical l=A. Upper set 
is for the observer perpendicular to the orbital plane and the lower set is for the observer 
in the orbital plane. Pictures show the stellar appearance at time intervals of 50 r g /c, with 
colour corresponding to the apparent temperature: gravitational redshift close to the black 
hole and Doppler shift of receding material stretch the observed frequency of photons (and 
therefore the observed temperature of the stellar surface) towards zero (red in colour code), 
while Doppler shift of approaching material increases the observed temperature (blue colour, 
corresponds to the value of twice the temperature in the system comoving with the star). 
Graphs show the apparent luminosity at different stages of the encounter (in units of the 
stellar luminosity before the encounter) and for different orbital angular momenta 1=0, 4, 5, 
7. (Full resolution images available at www.fmf.uni-lj.si/~ gomboc) 
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view from above the orbital plane 1^4 




time [tg/c] 

Fig. 5. — Star with R if =2r g and 1=4 during the encounter assuming no temperature change 
of the debris. Pictures show the stellar appearance at time intervals of 50 r g /c (except c and 
c', see graph) with colour corresponding to the apparent temperature according to the colour 
code: blue - temperature zero, white - 0.5 T c or higher. Inset graph shows the apparent 
luminosity at different stages of the encounter (in units of the initial luminosity far from 
the black hole) and for different orbital angular momenta /=0, 4, 5, 7. Upper set is for the 
observer perpendicular to the orbital plane and the lower set is for the observer in the orbital 
plane. (Full resolution images available at www.fmf.uni-lj.si/" gomboc) 
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lime [r e /c] time [r e /c] 

Fig. 6. — Effect of cooling of debris on the luminosity variations: solid line: no cooling, 
dotted: cooling by black body radiation, dashed lines: exponential cooling with decaying 
time 10 r g /c (long dash) and 1 r g /c (short dash). Results are for the star on 1=4 (left) 
and 1=5 orbit (right) as observed perpendicular to the orbital plane. (Full resolution images 
available at www.fmf.uni-lj.si/~ gomboc) 
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Fig. 7. — The effective eccentricity function of Roche penetration parameter. The 

lower three curves represent contributions due to the three excited modes (1,2,5). 
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